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Abstract. This work deals with the motion of a radially falling star in 
Schwarzschild geometry and correctly identifies radiation reaction terms by the 
perturbative method. The results are: i) identification of all terms up to first order 
in perturbations, second in trajectory deviation, and mixed terms including lowest 
order radiation reaction terms; ii) renormalisation of all divergent terms by the 
f Riemann and Hurwitz functions. The work implements a method previously 
identified by one of the authors and corrects some current misconceptions and 
results. 
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1. Introduction 

The burst of gravitational waves emitted during the capture of a compact star by a 
supermassive black hole is of main interest for space interferometry pQ, for testing 
general relativity in strong field and for investigating the physics of black holes. The 
detection of such sources requires the design of templates, which in turns require 
understanding of the complex orbital evolution during capture. These considerations 
are timely and strategical, especially if we refer to the supermassive black hole in the 
centre of our galaxy [2], as a source to be potentially detected in the next years. 
In the last years, perturbation methods have played an increasingly important role, 
due to their applicability to the last phases of lifetime of binaries constituted by 
comparable masses before merge [3]. The implication for ground interferometry 0] 
has not passed unnoticed [5] and considerable efforts and resources have been and are 
being employed. 

The understanding of orbital evolution requires understanding of radiation reaction 
effects. Indeed, radiation reaction, still partially outstanding problem in general 
relativity, like the two-body problem, is of most concern for gravitational waves 
detectors. Its influence is manifest in data analysis, where a phase mismatch of the 
templates with the signals may cause loss of detection. 

In this work, we analyse a radially falling star, m, captured by a massive black hole, M, 
by perturbative methods. The motion is studied in strong gravity, the perturbation 
being based on the m/M ratio. Radial fall is an idealisation of the capture scenario, 
but applicable to final plunging. Furthermore, most of the radiation, and thus reaction, 
occurs close to the horizon where inspiral has ceased. Nevertheless, we do not make 
specific reference to plunging as our analysis allows any distance for the point of 
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departure. 

We refrain from using energy balance and adiabatic hypothesis. The former is the 
imposition at start, rather than a rightful outcome, of the equality of the energy 
radiated with the energy loosed by the system. The latter can't be evoked since the 
particle immediately has to react to the radiation emitted, contrarily to inspiral motion 
where radiation reaction time scale is larger than the orbital period||. Furthermore, 
an analysis out of adiabacity implies that a straightforward linearisation of the 
phenomenon under study is not justified, and instead a careful screening of relevant 
terms must be adopted. Hence, neglection of higher order terms is to be avoided in 
absence of compelling reasons. 

We study the motion of a radially falling star with a perturbative approach using 
Moncrief gauge invariant formulation of the Regge- Wheeler [H] and Zerilli [§] 
equations. Incidentally, also second order perturbations formalism may be made gauge 
invariant |1U| . Under these assumptions, we shall not be concerned with gauge issues 
hereafter. 

Finally, second order perturbation analysis demands a consistent description 
at first order of the energy momentum tensor. The linearisation of general 
relativity implies that the particle motion in unperturbed Schwarzschild generates 
radiation; such radiation, including quadratic radiative terms of first order, plus the 
corrected motion, including radiation reaction, shall generate radiation at second 
order. Computing gravitational radiation to second order requires the knowledge 
of the trajectory of the falling mass on the first order metric (Schwarzschild plus 
perturbations). 

In the second section, it shall accordingly be developed a physical hierarchical scheme 
for each term constituting the geodesic, classifying terms that produce accelerations 
deviations from the unperturbed geodesic in Schwarzschild geometry in three types: 
the first depending upon the unperturbed metric evaluated on the perturbed trajectory; 
the second upon the perturbed metric evaluated on the unperturbed trajectory; 
the third upon the perturbed metric evaluated on the perturbed trajectory. In the 
third section, the issue of renormalisation shall be dealt for the infinite sum, on all 
multipoles, of finite terms that lead to divergencies. 

We adopt the geodesic concept in our approach, meaning that any motion is geodesic 
if the underlying metric is properly defined. 

Finally we are revising the work on the concept of self-force ^2] and correspondences 
between the two methods are subject of an undergoing investigation [13j . 

2. The metric, the perturbation scheme and the geodesic equation 

Perturbation method for analysis of radiation reaction has been previously proposed 
[TJ] - HE|. The metric is the sum of the Schwarzschild metric and the perturbations: 

4) >) -(i :,)~(t 

gtt = /(l - Ho) gtr = grt = -Hi g rr = - J (1 + H 2 ) 

|| Inspiral motion, especially around spinning black holes, at high eccentricity and on inclined orbits 
is characterised by three radiation time scales 
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9 U = y(l + #o) 5 tr = 9 rt = -Hi g rr = -/(l - H 2 ) 

r — 2M 

g»u = ffo, + K v = rf v -h» v / = 

r 

The order of the perturbation can be made explicit, i.e. h^\h^ 2 \ The 
position of the particle r e = r p + Ar p is given by the unperturbed trajectory in the 
unperturbed field r p and by several contributions, among which radiation reaction, 
given by the unperturbed and perturbed field, that generate a trajectory deviation 
Ar p . The field is developed in Taylor series around the real position of the particle: 
9fiu(r e ) — g^u{rp) + Ar p (dg^ v /dr) r . The geodesic is only dependent upon radial and 
time coordinates: 

^^^(i) 3 + (2 r.-r,)(i) 2 + (r.- 2 r : ;)(±Vr«, (D 



In absence of the weak field hypothesis, h is not limited in amplitude, but the 
following justifies that solely the terms in Tab. 1 are to be retained. We suppose: 

r^(l)i2 h {2) h (l) M 

^ < < 2 

T) T] T] r p 

The latter inequivalence is due to the twofold nature of the acceleration trajectory 
deviation: Ar p is the sum of two types of contributions. One is given by the 
Schwarzschild metric^, the other by the perturbations h (coupled with rj): 

Ar p = Ar p (r)) + Ar p (h) (3) 

The acceleration is dependent upon derivatives which are not necessarily 
small, especially in the last phase of the trajectory. In conclusion, the terms 
proportional to , Ar p , Ar p , derivatives and Ar^, Ar p , Ar p Af p are retained, 
while those to [/i^] 2 , and are neglected, as the second order terms in trajectory 
deviation when multiplied by first order perturbations. 

We write the geodesic equation in the following form: 

Af p = aiAr p + a 2 Afp + a 3 Arp + a 4 Afp + a 5 Ar p Afp + a 6 + a 7 Ar p + a$Ar p (4) 

The physical significance of the terms is essential + . The terms 01,2,3,4,5 arise 
from the pure Schwarzschild metric and they may be alternatively interpreted as 
representing the coefficients of the geodesic deviation of two particles separated on 
the radial axis by a Ar p distance. In the scenario of a single falling particle, they are 
the coefficients representing the unperturbed Schwarzschild metric influence calculated 
on the perturbed particle trajectory, i.e. the real position (0:1,3), velocity (02,4) or 
both (a^). The term is the lowest order containing the perturbations. It represents 
the perturbation influence calculated at the position and velocity of the particle in the 
unperturbed trajectory. Thus it does not represent radiation reaction, although it 
contributes to. The 07 term represents the perturbation influence calculated on the 
real position of the particle in the perturbed trajectory, h^Ar p . Finally, the as term 
represents the perturbation influence calculated on the real velocity of the particle in 

% Indeed, the first five terms a\— 5 provide an acceleration f p depending solely upon 77, the 
unperturbed metric 

+ The ai term correspond to the A term of |17l . |18| . apart of an error in the quoted publications 
(see appendix). The 02 term correspond to the B term, ctg correspond to C of 1171 . 1181 . 
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the perturbed trajectory h^ v Ar p . 

It is thus a development at first order in perturbations and second order in trajectory 
deviation or mixed terms. Tab. 1 lists all a terms. We emphasize that: 

• The terms 03,4,5,7,8 were previously [T7| . [T5] neglected and the remaining terms 
a i,2,6 do not really represent radiation reaction; ii) the terms 03,4,5 are of second 
order in trajectory deviation. They also do not represent radiation reaction but 
are kept for mathematical consistency with the hypothesis The terms 017,8 
represent the lowest order radiation reaction terms. 

• Solely terms in second order of perturbations are excluded. Such choice is also 
compliant to the nature of the energy-momentum tensor of the second order 
equations. The tensor T^ v , that generates second order perturbations, must entail 
a geodesic stemmed from eq.(0J. The linearisation of general relativity induces 
a stepped-up approach, in where radiation and motion assume sequentially 
the lead. Motion generates radiation that generates reaction, which in turns 
generates radiation of second order (to be dealt with squared first order radiation 
terms). But radiation of second order requires generation by radiation reaction 
at first (lowest) order, i.e. the terms 0:7,8- 

• Radiation reaction effects, supposedly, are larger near the horizon, where most 
of the radiation is emitted. The leading terms of a§ near the horizon are of 
the order h rr r) rr ^ r r p + rf t h u ,t'i'p + h rr rj t t,r, whereas in 07,8, are of the order 
h TTr lrr,rr r p + rflhtt,t r p an d h rr rj rrr r p + r\ u h ttt to be multiplied by Ar p and 
Ar p , respectively. The Schwarzschild metric tensor components and derivatives, 
above quoted, tend to infinity and its powers, while r p to zero, near the horizon. 
Estimate of the leading terms in 017.8 justifies their inclusion. 

• Second order perturbations equations, are yet unsolved in presence of a source 
term and the vacuum solutions are adequate for dealing the "close limit" but 
not the particle motion. The unavailability of the second order perturbation 
solutions suggest to consider the problem pragmatically using all remaining known 
quantities, as shown in eq.Q. 



3. Black hole polar perturbations equation 

Zcrilli 9 found the equation for polar perturbations and studied the emitted radiation 
adding a source term, a freely falling test mass m into the black M. The equation 
is written in terms of the wavefunction ipi for each 1-pole component, the tortoise 
coordinate r* , the polar potential Vi(r), the 2 z -pole source component Si(r,t): 

d 2 M r ,t) d 2 ^, t) „ / r \ 

—^2 ^2 Vi(r)Mr,t) =Si(r,t) r = r + 2Mln ^— - I J (5) 

w - (1 - g) 2A2(A + iy %ZmT m ^ + 18MS H™><°) 

1 - ^M}j U4y/(2l + 1)tt 
Sl = (A + l)(Ar + 3M) X 

It shall be the upcoming numerical simulation to verify the relative weight of these terms 1191 . 
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2MV 



S'[r-r p (t)] 



x , M 

X+l 

r 



6Mr 
Ar + 3M 



*[r-r p (t)] 



where r p (£), geodesic in unperturbed Schwarzschild metric, is the inverse of: 

-l" 



-AM 



/ r \ V 2 4M / r \ 



3/2 



\2M) 



3 \2M J 



2M In 



2M 



- 1 



The perturbations around the particle are (Regge- Wheeler gauge H = H 2 



r 
2M 

i 



1 



9M 3 + 9AM 2 r + 3 A 2 Mr 2 + A 2 (A + l)r 



r 2 (Ar + 3M) 2 

3M 2 - AMr + Ar 2 
r(Ar + 3M) 



■0 + 



V», r + (r - 2M)^>, 



, , 3Af 2 + 3AM r - Ar 2 , 
ffl = (r-2M)(Ar + 3M)^- 4 

The unperturbed velocity is given by (ro is the test mass position at start): 



(7) 



(8) 



2M 



2M\ 
ro J 



1/2 



2M 
ro 



-1/2 



(9) 
(10) 

(11) 



4. Renormalisation 

In this section we reconstruct the renormalisation for ae and apply, for the first time, 
the £ function to a^g. Indeed, the infinite sum over the finite multipole components 
contributions leads to the problem of dealing infinities in the results. For ever larger I 
the metric perturbations tend to an asymptotic behaviour. In other words, the curves 
representing each metric perturbation component for each 1, accumulate over the 
I — > oo curve. Thus the subtraction from each mode of the I — > oo leads to a convergent 
series. We extend the application of the Riemann £ function for renormalisation 
|T7) - ^B] to all pertinent terms of the geodesic of Tab. 1. Instead, mode-sum 
renormalisation is planned in the near future. For L = I + 0.5, the wavefunction 
and its derivatives assume the following forms at large L or I [T^], 120111 when averaged 
around the particle at r p : 



1p,r 



6V27rm(r - 2M) 
r {r p -2M) 



L 



-2.5 



4V27rm(r - 2M) 



6V27rmVn) - 2Mr 



ror p 



r (r p -2Mf 
P L- 25 
4\^2nmy/ro 



5(r - 2M) 



9M 



2r Q 



- 6 



4V27TTO(r - 2M) 
r (r p -2M) 2 

£-0.5 



4V27rm v / ro - 2M rj 



2Mr'„ 



ror p (r p - 2M) 5 



?torp(T- p - 2M) 
5(r - 2M) 



-L 



-0.5 



2r n 



9M 



L 



-0.5 



J The derivation of such expressions, quoted from a paper in preparation by Barack and Lousto, as 
referred by 1181 and 1211 . is yet unpublished. 
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Using the above equations, eas. (|9H0|l and recasting H l t as function of 
tjj, ip r , ip rr , ip rrri the term for large L or I around the particle is: 
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£< 

1=0 



a l e = a a & L° + a\L- 2 + alL-* + 0{L-*) (12) 



The term a^L needsff renormalisation 22 . The Riemann ( function and 
its generalisation, the Hurwitz £ function [23j, are defined by: 



i=o 



1=1 

where in our case a = 0.5. Thus: 

oo 

C(s,0.5) = ^0 + 0.5r s = 2 s 



i=o 



^2(2i + iy 



1=0 



Due to the imparity of the term in braces, ea. H14(l is rewritten as: 



C(s,0.5) = 2 s {((s)- 



.1=0 



= 2 s (1 - 2- s ) C(s) = (2 s - l)C(s) 



(13) 



(14) 



(15) 



Some special values of the Hurwitz functions are: 



1 



C(2,0.5) = -^ 2 



C(4,0.5) = -tt 4 (16) 



C(-2,0.5) = C(0,0.5)=0 

The latter values when applied to ea. (|12f) . give: 

oo oc oo 

a 6 = a a 6 + 0-5)° + 4 + °- 5 )" 2 + < £(' + °- 5 )" 4 + + °- 5 )" 6 ] 



1=0 



1 = 



agC(0,0.5)+^C(2,0.5)+a^C(4,0.5) + [0(;+0.5)- 6 ] = I 7r 2 a fc+jU 4 ag+[0(/+0.5)- 6 ](17) 

For the renormalisation of a^^ terms, the expressions: ip jtip jttip jtrip ,trrrip ,rrrr 
are deducted |22 operating on the averaged wavefunctions and derivatives, and the 
homogeneous wave equation. The latter is recast as ^31 - |15| : 

1 d 2 Mr,t) cP^M) , p - 1 dil>i{r,t) 



dr 2 



dt 2 



rp 2 dr 



Vi(r)il>i(r,t) = (18) 



where p = dr* /dr. Deriving sequentially ea. i|18|) . we get the tj) needed derivatives. 
The latter are evaluated for L — > oo and when inserted in the ay s terms, result into: 



a 7 = £ 4 



1=0 



£< 

1=0 



a\ = a a 7 L 2 + a b 7 L° + a c 7 L~ 2 + afL~ 4 + 0(iT 6 ) 



4 = a%L° + a\L- 2 + o^L^ + 0(iT 6 ) 



(19) 
(20) 



Renormalisation of eqs. i|19l2U|> leads to: 



a? = \^a 7 + \^a d 7 + [0(1 + 0.5)- 6 ] a 8 = \v 2 a\ + -7r 4 a| + [0(/ + 0.5)- 6 ] (21) 
2 o 2 o 

ft The term a%L° differs from b of eq.(13) in [TS| which is again different from the value given in 1201 . 
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5. Conclusions 

We have obtained the following results: i) calculation and determination of all terms 
up to first order in perturbations, second in trajectory deviation, mixed term of second 
order including lowest order radiation reaction terms, all contributing to the trajectory 
of a radially falling test mass in Schwarzschild geometry; ii) renormalisation of all 
divergent terms, including new ones, stemmed from the infinite sum of finite angular 
momentum dependent components by the £ Riemann and Hurwitz functions; iii) 
correction and improvements of previously published results. 
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Appendix on the geodesic equation 

Three perturbations schemes concur: 

r e = r p + Ar p g pv (r e ) = g pv {r p ) + Ar p (^-jjfrj 9^ = f» + hpv (22) 

Hence, the development of eq. Q leads to a lengthy computation 
/ dr\ 3 1 

r ' r I dt J ~ 2 [ gU ( 29tr > r ~ 9rr ^ + 9 tr 9rr, r ] (f| + 3r 2 p Ar p + 3r p Af;) 
i [(77" + h u + n u r Ar p + h u r Ar p ) {2h tr , r + 2h tr , rr Ar p - h rr>t - h rrttr Ar p ) + 



(h tr + h^Arp) (r)rr,r + h rr ,r + r] rr ^ r Ar p + h rr , rr Ar p )\ (r p + 3f 2 Ar p + 3r p Ar 2 ) (23) 
2r* r (^) = [g tt g tt ,r + g tr grr,t] (r 2 p + 2f p Af p + Ar 2 p ) ~ 

[(77" + h u + ri u r Ar p + h u r Ar p ) (r) tt , r + h tt , r + Vtt,rrAr p + h tt , rr Ar p ) + 
(h tr + h%Ar p ) (h rr , t + h„, tr Ar p )] (r 2 p + 2f p Ar p + Ar 2 p ) (24) 
/ dr\ 2 1 

— [(?7 rr + h rr + rfj Ar p + h r ^Ar p ) [r} rTjr + h rT%r + r} rrj „Ar p + h rr ^ rr Ar p ) + 

(h rt + h rt r Ar p ) [2h tr ,r + 2h tr . rr Ar p - h rr<t - /i rr , tr Ar p )] [r 2 p + 2r p Ar p + Ar 2 p ) (25) 
( dr \ 1 

tt [ dt J = 2 t 5 " 5 "^ + 5 (2ffrM ~ 5 "' r) ] ^ p + ^ fp) ~ 
i [(??" + h** + r;"Ar p + h u r Ar p ) (h tt ,t + h tt , rr Ar p ) + (26) 

(h tr + h^Atp) (2h rt ,t + 2h rt: trAr p - T]u,r - htt,r - T)tt,rrAr p - h u ,rrAr p )] {f p + Af p ) 

-2T r tr = - [g rr g rr , t + g rt g tt ,r] (r p + Ar p ) ~ 



- [{rf r + h rr + r)JAr p + h r rAr p ) (h rr , t + h rrttr Ar p ) + 

(h rt + h ri r Ar p ) (r) tt ,r + h tt , r + r) tt ,r r Ar p + h tttrr Ar p )] (r p + Af p ) (27) 

-Hi = -\ [9 rr (2ff rt ,t - gtur) + g rt gu,t] ^ 

~\ [(^ + h " + ^"Arp + h r JAr p ) (2h rt , t + 2h rtjtr Ar p - r) tt ,r ~ h t t,r - iltt, rr Ar p ~ h u ,r r Ar p ) 
+ (h rt + h«Ar p ) (h tt , t + h tt , rr Ar p )] (28) 
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The development of eas. (|23l - l28[) produces the coefficients vertically listed in Tab. 
1. The lines, hierarchically developed in horizontal rows, show the a,; (i = 1,..8) terms 
of eq.(@J. The e.g. a.\ term is given by: 



1 



1 



a l = tfrVtt^rp + ■q Ll 'qtt,rrrp - r^'V,^ ~ ^"^"'^P + ^rVtt,r + ^V^Vt^rr (29) 

and results into ea. (|31 \ , It's obvious that a similar development takes place for 
all a terms. But for the numerical estimate of aej,$ terms, simulation of the radial 
fall in time domain is required 19 . The simulation is to provide wavefunctions and 
thus metric perturbations. Here below all a terms are shown (Regge- Wheeler gauge 



a 



Oi\ = 



M 
r 

2M 



2M 



3M 



3M 



-2M p 
3(r — M) 



1 



r (r- 2M) 



2 V 



a 2 



a 3 



Q'4 



a 5 



6M 



r(r- 2M) p 



AM 2 



AM 



(r - 2M) 



3 P 



3M 



r{r - 2M) 

AM 
~r(r - 2M) 

1 



M 



«(> 



+ 



r - 2M 
r-2M 



r(r - 2M) 

„2 



2(r - 2M) 



H - 



2 

- + 

r 

M 



1 



r-2M 



- 2M 
Hi - rH[ 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



H' f 2 p -3[ -H - '—Hi ] r, 



1 • 



M 



2M 



a 7 



2M 



2Mr 



5 H" r 2 - - 
2 p 2 



(r - 2M) 2 
AM 



Hn--- 



H' - Hi 



1 



(36) 



H' 



2r-2M 
2M A 



:H'- 



2M 



(r - 2M)' 



M 



2M 



r 3 
p 



AM(r — 3M) AM(r-2M) l(r-2Mf 

~3 ^0 T2 ^0 - ^ ; #0 



(37) 



2M ■ 

H x + (r - 2M)H[ 



«8 



2M 



1 i 

2~2M' 



-H n - 



M 



r - 2M 



-Hi - rH{ 



r 2 
p 



3H' r p 



-Hq H ttHi 



(38) 
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Table 1. Classification of geodesic terms. 



**-(£)" --m)' n.(S) -n, 







1 rr -2 
2^7 r]rr y rrp 




\rf r r\tt,r 




Unperturbed Schwarzschild ao term 








rfrVtt, r fl 


1 rr -2 
2^/) r rjrr,rrp 








rj U ntt,rrfl 


1 rr -2 
2~V r/rr,rrrp 

qi term 




\rf r r)tt,rr 




2rf t r)tt,rr p 


V r/r^rp 

Q2 term 








r) u r r\tt,rrr 2 v 


2^?,r rjrr.rrrp 

Q3 term 




\rf^r]tt,rr 




rj n t t, r 


1 ^rr 

2 w r\rr,r 

ol^ term 










V,r r/rr,rrp 








! ir\ u r\u,rvrp 


*H rjrr,rrrp 

a 5 term 






V U htr,rfp 


-h u r)tt, r rl 


\h rr rirr,rrl ^V Uh tt,tfp 


?^ hrr,t^p 


-rf r h r t,t 




V U htt, r rl 


-\ri rr hrr,rfp \h U ' Vtt , r f p 


h rt ritt,rf p 


-\h rr nu,r 


-lh tr n f 3 

2 't lfrr,rl p 




«6 term 




1 rr l 
2» /ltt,r 




-h u r r\tt,Tr 2 v 


1 7 rr ■ 2 1 tt l. 

2 > r ^]rr,r^p 2^> r **j*^"p 


^?,r ^rr^t^p 


—rfjh r t,t 
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